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Abstract. The standard analysis of the CMB data assumes that the distance to 
the last scattering surface can be calculated using the distance-redshift relation as in 
the Friedmann model. However, in the inhomogeneous universe, even if (dp) = 0, 
the distance relation is not the same as in the unperturbed universe. This can be of 
serious consequences as a change of distance affects the mapping of CMB temperature 
fluctuations into the angular power spectrum C'i . In addition, if the change of distance 
is relatively uniform no new temperature fluctuations are generated. It is therefore 
a different effect than the lensing or ISW effects which introduce additional CMB 
anisotropics. This paper shows that the accuracy of the CMB analysis can be impaired 
by the accuracy of calculation of the distance within the cosmological models. Since 
this effect has not been fully explored before, to test how the inhomogeneities affect the 
distance-redshift relation, several methods are examined: the Dyer-Roeder relation, 
lensing approximation, and non-linear Swiss-Cheese model. In all cases, the distance 
to the last scattering surface is different than when homogeneity is assumed. The 
difference can be as low as 1% and as high as 80%. Excluding extreme cases, the 
distance changes by about 20 - 30%. Since the distance to the last scattering surface 
is set by the position of the CMB peaks, in order to have a good fit, the distance needs 
to be adjusted. After correcting the distance, the cosmological parameters change. 
Therefore, a not properly estimated distance to the last scattering surface can be a 
major source of systematics. This paper shows that if inhomogeneities are taken into 
account when calculating the distance then models with positive spatial curvature and 
with S1a ~ 0.8 — 0.9 are preferred. The ACDM model (i.e. a flat Friedmann solution 
with the cosmological constant), in most cases, is at odds with the current data. 



1. Introduction 

It has been said that we entered the era of precision cosmology. This is mainly due to 
observations of the cosmic microwave background radiation (CMB). As the CMB power 
spectrum is very sensitive to cosmological parameters it provides very tight constraints. 
Currently the errors are at the level of a few percent [lj. The CMB power spectrum 
is shaped by several processes that can be divided into two groups. The first group 
involves processes that occurred before and during the generation of the CMB. This 
part is well understood as the Universe at that times was very close to homogeneous and 
contribution from the spatial curvature and dark energy was negligible. The constraints 
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on cosmological parameters coming just from this group of processes ware recently 
presented and discussed in [2J. The second group of processes takes into account what 
happens with photons between the last scattering surface and the observer. This includes 
the distance to the last scattering surface, the ISW effect, the impact of reionization on 
the CMB, the Sunyaev-Zel'dovich and the lensing effects. 

The analysis of this second group of effects is a subject to systematics, as the 
geometry of the late time Universe does not have to be close to homogeneous and 
isotropic Robertson- Walker geometry (these processes are most often studied within the 
framework of linear perturbations around the homogeneous Friedmann models). On 
small scales (say I > 60, i.e. 9 < 3°), expect for the distance, these processes should 
not significantly modify the overall shape of the CMB power spectrum. The distance 
to the last scattering surface affects the mapping of CMB temperature fluctuations into 
the angular power spectrum C\. This effect can be modelled by introducing the shift 
parameter 71, which is proportional to the distance to the last scattering surface [H 0]. 
Therefore, there is a number of papers where the constraints from the CMB are just 
limited to the constraints on the shift parameter. This is usually done when testing 
different models of dark energy (5} |6], alternative cosmological models [H |8], or Gpc- 
scale inhomogeneous models |2l|TDlinilI2llI3llIlllI5llISlinilIB]. As discussed in [HI [20] 
this type of analysis can be improved but taking into account both the scale of the sound 
horizon at the last scattering and TZ. 

This paper studies the effect of the distance but does not consider any alternative 
cosmologies. Here we just focus on the effect of small-scale inhomogeneities on the 
distance to the last scattering. The presence of inhomogeneities is know to affects 
the distance-redshift relation [2T1 [22 j. Some studies claim that the effect is large 
[231 121 [251 HE], others that it is small [271 12H1 [291 EH [31] • Even if the change of the 
distance just at the level of a few percent [32] [33j [3H [35] this analysis is still important 
as a change of the distance by a few percent leads to a similar change of the inferred 
values of cosmological parameters. 

One may think that the effect of inhomogeneities is taken into account in the 
standard analysis of the CMB via the ISW or lensing effects. However, in both cases 
(in the standard approach) this is done by using the matter power spectrum. Thus, this 
analysis is insensitive to the change of the mean, i.e. the uniform change of the distance 
to the last scattering surface. For example the lensing analysis deals with the change 
of the distance but is only sensitive to the change of the variance, not the mean (cf. 
[361 137]). If the change of distance is relatively uniform (i.e. with a negligible variance) 
then no additional temperature anisotropics are generated. Since the change of distance 
affects the mapping of the physical position of the peaks to the peaks in the angular 
power spectrum Ci this effect needs to be considered also at the level of the mean. 

Therefore, if inhomogeneities are not taken into account when calculating the 
distance, then they may become a major source of systematics. With the increasing 
precision of CMB experiments, in particular Planck, a proper understanding of 
systematics is important. Without taking into account all systematics the precision 
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cosmology will not be an accurate cosmology. Therefore, this paper aims to study how 
the presence of inhomogeneities affects the distance-redshift relation, in particular the 
distance to the last scattering surface. Several models are considered and it is shown 
that the change of distance due to inhomogeneities is not negligible. 

The structure of this paper is as follows: Sec. |2] discusses different methods for the 
distance calculation, Sec. [3] applies these methods to calculate the distance to the last 
scattering surface, Sec. H] explores the implication of differences in the distance, and 
Sec. [5] concludes the results. 



2. Distance and distance corrections 



As shown by Sachs |21| . the equation for the angular diameter distance D A is 

^ = -(M 2 + ^/^)A4, (1) 

where s is a parameter along the geodesic, k a is a null vector that is tangent to the 
geodesic, a is the shear of the light bundle, R a p is the Ricci tensor and R a ^k a k^ = 
KT a pk a kP. In the comoving and synchronous coordinates, for pressure-less matter, 
Rap^k 13 = npk k . As seen, the distance depends on density fluctuations, and secondly, 
even if the perturbations vanish after averaging (i.e. (p) = (po + 5p) = po, where po 
is the background density) they do modify the distance relation and the final result 
deviates from the homogeneous solution. Equation (CQ) is the differential equation so it 
is sensitive to any change of initial conditions. Therefore, it is sensitive to 5pD A as this 
term acts as a potential (for a detailed discussion see |3l]). As the parameter s is not 
directly measured (unless the geodesic parametrization is chosen in a particular way), it 
is convenient to express it in terms of a redshift. The redshift is given by the following 
relation 138 



l + g = p^Ql } (2) 

where u a is the 4-velocity, and subscripts e and o refer to emission and observation 
instants respectively. In the comoving and synchronous coordinates u a = 5$ and 

where k° = dt/ds. Thus, for a given model and using the above relation the parameter 
s can be replaced with z. 



2.1. Homogeneous universe 

When homogeneity is assumed, i.e. a = = 5p, and k° = 1/a (where a(t) is the scale 
factor) then dz/ds = —(1 + z) 2 H ^Q m (l + z) 3 + Q k (l + z) 2 + Q A and ([!]) reduces to 



D A ( Z ) = — —= sin I I' dz' , = ) , (3) 
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where H is the Hubble constant, Qj, = —k/(HQao), Q m = (87rGp ) / '(3-ffg), and 

q a = a/(3h 2 ). 

2.2. Dyer-Roeder approximation 

The Dyer-Roeder approach also assumes homogeneity but takes into account that light 
propagates through vacuum. Therefore, Q m which photons 'feel' is different than true 
Q m . This is modelled by a constant parameter a (of value between and 1) that 
multiplies Q m . In this case ([!]) reduces to [391 HQj I4"T] 



d 2 D A ( 1 dH 2 \ dD A 3H 2 M 

-d^ + {H^ + TT- z )^ z - + 2 S^ 1 + z)Da = °' (4) 

where H(z) = H ^yil m (l + z) 3 + 0^(1 + z) 2 + O^. In order to include the evolution of 
density fluctuations along the line of sight, Ref. |35| suggests the following modification 
of the Dyer-Roeder equation 

^ U^JL^) ^H-fifU^) (l + Da — 0, (5) 



dz 2 \H dz 1 + zJ dz 2H 2 my '\ {l + zfl\ 

where {S)\b is the mean of the present-day density fluctuations along the line of sight. 
Initial conditions needed to solve (@J and (jSJ) are: Da(0) = and dD A /dz = l/H . 

2. 3. Linear perturbations and the lensing approximation 

Within the linear perturbative scheme the distance is 

D A (z) = D A (l + A D ), (6) 

where D A the distance in the homogeneous universe and is given by ([3]). The most 
general form for is presented in [4^| |4"3| [HJ US]. Neglecting the contribution from 
the observer's and source's motion, and taking the leading term, can be expressed 
as 

Xe 

A D = - / d X ^^ X V 2 0(x), (7) 

J Xe 


where x is the comoving coordinate d\ = dz/H(z), <fi is the gravitational potential 
that relates to the density perturbations 5p via the Poisson equation V 2 = ^£-a 2 5p. 
The above expression is equivalent to the convergence in the lensing approximation and 
is known as the Born approximation. As seen, propagation through voids (dp < 0) 
increases the distance while propagation through overdense regions (8p > 0) decreases 
the distance. 

In order to solve ([7j) one needs to know exactly what the density fluctuations along 
the line of sight are. However, within the weak lensing analysis one does not solve 
directly. Instead the variance is calculated, thus instead of Sp one has (Sp) 2 which can be 
expressed in terms of the matter power spectrum. In this paper we are interested in the 
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overall (uniform) change of the distance. Thus, we will solve (J7|) directly, this requires 
knowledge of density fluctuations along the line of sight. The algorithm is described 
below, but see also |l6j H7] for an alternative approach. 

3. Distance to the last scattering surface 

To test what is the effect of inhomogeneities on the distance to the last scattering surface 
and the accuracy of the CMB analysis, let us consider 7 different ways of calculating 
the distance to the last scattering surface 

(i) The ACDM model: 

Distance to the last scattering surface is calculated from ([3]). 

(ii) Dyer-Roeder approximation (model DR): 

Based on the cosmological observations, Yu et al. [48j suggest that a = 0.93. 
Therefore, the distance to the last scattering surface is calculated using (0J with 
a = 0.93. 

(iii) Modified Dyer-Roeder relation (model mDR): 

The modified version of the Dyer-Roeder equation (j3J) is solved with (S)id = —0.3. 

(iv) Lensing approximation (model Len 1): 

In order to solve ([7]), one needs to know V 2 </> along the line of sight. Using the 
Poisson equation the gravitational potential is related to density fluctuations. The 
density fluctuations along the line of slight are generated from the log-normal PDF 
- the algorithm is discussed in details in Appendix A 

(v) Lensing approximation (model Len 2): 

In model Len 1 the mean of density fluctuations along the line of sight is negligible. 
However, as pointed out in [35] the mean of density fluctuations along the line 
of sight does not have to vanish, as (5)sd = (3D - volume average) does not 
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necessarily imply (5) id = (ID - along the line of sight). In model Len 2, (5) id ^ 



although (6)30 — 0. The algorithm is discussed in details in Appendix B 
(vi) Lensing approximation (model Len 3): 

Model Len 3 uses a more flexible (than Len 1 & 2) method of generating density 
fluctuations along the line of sight. The algorithm is presented in details in 



Appendix C 



(vii) Fully non-linear model (model SC): 

The above approaches were based either on the assumption of homogeneity (the 
ACDM and Dyer-Roeder relation) or on the linear approximation (models Len 
1-3). In fact, the lensing approximation is not even fully linear, as it takes into 
account just the leading linear term. Also, when generating fluctuations along 
the line of sight, the matching conditions were not taken into account, i.e. the 
structures were not properly matched. In addition, the light rays were not properly 
match as well, i.e. when light ray exits one structure and enters another one, the 
null geodesies need to be properly matched. Therefore, a fully non-linear model is 
considered here, where the distance is calculated directly from ([T|) and all matching 
conditions are handled properly. The model is based on the Swiss-Cheese model. 
Each inhomogeneous patch is of spherical symmetry and is described using the 
Lemaitre-Tolman (LT) model [1H1 [50] . Thus, density within each inhomogeneous 
patch changes continuously, unlike in the lensing approximation, where within each 
generated structures, matter density is constant. Also the evolution of matter is 
calculated using the exact method, not the linear approximation. The properties 
of the model are chosen so that the mean of density fluctuations along the line of 
sight is zero, i.e. (5) id ~ 0. Thus, this model allows us to test how large non-linear 
contributions are, in contrast to (5) id 7^ contributions. The detailed description 
of this model is presented in Appendix D 



In this section it will be assumed that the large-scale universe is exactly the ACDM, 
so Hq = 71.4 km s _1 Mpc -1 , Q m = 0.262, Q\ = 0.738. A sample of density fluctuations 
along the line of sight is presented in Figure [TJ The results, in the form of Da(z), 
are presented in Figure [2j The distance to the last scatting surface is presented in the 
inset, and in the form of the shift parameter 1Z in the 2nd column of Table [TJ As seen, 
the distance to the last scattering surface depends on the choice of the method. Even 
within the most optimistic case, model Len 1, the change is of order of 1%. However, the 
non-linear contributions (even if (5) id ~ 0) gives ~ 30% change of the distance. When 
(3) id 7^ 0, even within the lensing approximation the change is of order of 20 — 30%. 

Another important fact is that the variance within each method is small. Models 
Len 1-3 allow for random (within the limits of the method) density fluctuations along 
the line of sight. Based on 100,000 runs (each with a different distribution of density 
fluctuations along the line of sight) the results are Len 1: = —0.0071 ± 0.0198, Len 
2: A D = 0.2033 ± 0.0064, Len 3: A D = 0.2729 ± 0.0071. This shows that the change 
of distance is relatively uniform. Below, when referring to models Len 1-3 it will be 
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Figure 2. The angular diameter distance within 7 considered models. The inset 
presents Da in the vicinity of the last scattering surface. 



referred to the mean of Ad |j. 

In the next section we will study the implications of the change of the distance 
to the last scattering surface and the CMB constraints when this change is taken into 
account. 

4. CMB analysis with the adjusted distance-redshift relation 

4-1. Fitting the distance and the accuracy of CMB analysis 

In this section we assume that the physical densities of baryonic and cold dark matter, 
and primordial power spectrum are the same as in the standard cosmological model, 
and only implications of the change of the distance to the last scattering surface are 
going to be examined. 

The change of the distance-redshift relation leads to a different mapping of the CMB 
temperature fluctuations into the angular power spectrum C\. Thus, if the distance- 
redshift relation changes, then the standard analysis of the CMB needs to be adjusted. 
The way to do this is following: first the CMB power spectrum is calculated using 
standard codes (for example CAMB [51J) then the output requires following scaling 

£^£(1 + A D ) and C e C t (l + A D y 2 . (8) 

Thus, the £ axis changes while the amplitude £[£ + 1)C^ remains relatively unchanged 
(for a detailed discussion see [T2l El HH]). To present how it works let us consider the 
following example: let us change the value of A in such a way so that the distance to the 
last scattering instant, for studied models, is the same as in the standard cosmological 

X It should be noted that the distance corrections depends on cosmological parameters of the 
background model, for example Ad for model Len 1 changes to —0.0064 and —0.0069 when parameters 
are as in Table [1] and Table [2] respectively. 
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model. Q\ needed to fit the distance, for studied models, is presented in the 3rd column 
of Table [TJ Apart from model Len 1, Q\ is much larger than when homogeneity is 
assumed. If the power spectrum was calculated for a homogeneous Friedmann model 
(i.e. with the distance-redshift relation the same as in the Friedmann model) the result 
would look like the one presented in Figure [3j However, after applying scaling (jBJ) the 
angular power spectrum is as in Figure HI As seen the CMB power spectrum is very 
similar except for low-£ where the ISW effect is different. The low-£ power spectrum is 
presented in the inset, where the error-bars are enlarged to include the cosmic variance 
AQ/Q = ±(2/(2£+l)) 1 /2. 

Another example: let us correct the distance by changing the cosmological 
parameters but keeping the spatial curvature zero, i.e. Q m + = 1. Thus, the change 
of fl m automatically implies the change of Q\. Also to keep the shape of the CMB 
power spectrum unchanged, Vl m h 2 = (Qb + Q c )h 2 must not change, so the change of fl m 
automatically implies the change of H . The change of cosmological parameters needed 
to fit the distance is presented in Table [3 As seen all models except for model Len 1 
can be ruled out because of too high amplitude of H . Even within model Len 1, the 
change of parameters is large compared to the current precision - for example in model 
Len 1, Q m is around 0.28 not 0.262 so the difference is around 7%, which is twice as 
large as the error estimated by the WMAP team [lj. For completeness the CMB power 
spectrum for Len 1 and the ACDM model are presented in Figure [51 

The above examples show why it is important to have an accurate method for 
calculating the distance. As the universe is inhomogeneous, the distance-redshift relation 
is different than the distance-redshift relation in the Friedmann model. As shown, even 
a small change of the distance can have a large impact on the accuracy of the CMB 
analysis. As seen from Tabled] a change of the distance by 1% leads to the 0.5% change 
of the cosmological constant. Although this change is small, in the era of precision 
cosmology this may be important. Especially that as shown in Table |2] if one puts by 
hand the assumption of spatial flatness, then systematics are much larger. 

4-2. Spatial curvature of the Universe 

Let us now focus only on models Len 1, Len 2, and SC. We exclude the Dyer-Roeder 
relation as this is only a rough approximation (as the shear accumulates along geodesies 
this approximation is expected to be inaccurate for large distances). Also as model 
Len 3 is rather of arbitrary matter fluctuations (it was considered only to check how 
arbitrariness of matter fluctuations changes the results) it is excluded from the further 
and more detailed analysis. 

In this section we will only focus on Len 1, Len 2, and SC models. Len 1 model 
is a conservative approach - it is based on the linear approach with random density 
fluctuations along the line of sight. Len 2 model is also based on linear approximation, 
but the mean of density fluctuations along the line of sight is non-zero (as discussed 
in [35] if structures in the universe are organized, then density fluctuations are no 
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Table 1. The effect of inhomogeneities on the distance to the last scattering surface 
- 2nd column (TZ = Da(1 + z)y / fl m Ho/c); the accuracy of estimation of A - 3rd and 
4th column. In all models H Q = 71.4 km s' 1 Mpc' 1 , il c h 2 = 0.1107, VL b h 2 = 0.0227. 



Model the shift parameter 7Z Vt A needed to 5\ = |A^a|/^a 
(if the background model fit the CMB 
is the ACDM model) 



ACDM 


1.73 


0.738 


0.0 


Len 1 


1.71 


0.734 


0.005 


Len 2 


2.08 


0.797 


0.08 


Len 3 


2.19 


0.802 


0.09 


SC 


2.31 


0.919 


0.24 


DR 


3.19 


1.071 


0.45 


mDR 


2.04 


0.823 


0.12 




Table 2. The effect of inhomogeneities on the accuracy of inferred cosmological 
parameters. In all models il c h 2 = 0.1107, ft b h 2 = 0.0227, and Vt m + Oa = 1. 



Model 






n A 


5 A = |A^ A |/n A 


h 


S h = \Ah\/h 


ACDM 


0.262 





0.738 





0.714 





Len 1 


0.281 


0.07 


0.719 


0.03 


0.690 


0.03 


Len 2 


0.050 


0.81 


0.95 


0.29 


1.64 


1.30 


Len 3 


0.034 


0.87 


0.966 


0.31 


1.982 


1.78 


SC 


0.024 


0.91 


0.976 


0.32 


2.344 


2.28 


DR 


0.006 


0.98 


0.994 


0.35 


4.568 


5.40 


mDR 


0.095 


0.64 


0.905 


0.23 


1.19 


0.66 
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Figure 4. The CMB power spectrum calculated with CAMB [51| and corrected for the 
distance (i.e. with the scaling (jSJ employed). The inset presents the power spectrum 
for \ow-£ with errors enlarged to include the cosmic variance. 
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Figure 5. The CMB power spectrum calculated with CAMB [51] and corrected for the 
distance (i.e. with the scaling (j8|) employed). The inset presents the power spectrum 
for \ow-£ with errors enlarged to include the cosmic variance. 



longer purely random, hence even though the 3D average vanish, the average of density 
fluctuations over the the line of sight is not zero). SC model is a fully non-linear model 
with vanishing density fluctuations along the line of sight. 

The results of the CMB analysis within these 3 models is presented in Fig. [6j Now 
we allow Q m and to vary freely (however, we keep fl^h 2 and Q c h 2 unchanged, so the 
change of fl m implies the change of H ). As seen the constraints obtained within Len 
1 model are very similar to the constraints obtained within the standard cosmological 
model. However, model Len 1 does not include non-linear corrections. If these are taken 
into account or if only the mean of density fluctuations along the line of sight does not 
vanish then spatial curvature cannot be zero. As seen the presence of small-scale non- 
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Figure 6. The CMB constraints after taking into account the effect of inhomogeneities 
on the distance to the surface of last scattering, for Len 1, Len 2, and SC models. 
The color regions are the 95% confidence interval, the yellow lines indicate the 68% 
confidence interval. 



linear structures changes the distance-redshift relation and now the CMB data combined 
with local value of H no longer favors spatially flat models. This shows that the flat 
Friedmann solution is not preferred by the data, and when inhomogeneities are included 
into the distance-redshift relation then k > models with larger amount of dark energy 
are favored. 

Let us discuss for a moment the possibility that in fact Q\ ~ 0.8. What are the 
implications? In this case, the age of the universe increases to approximately 14 Gyr. 
This change is acceptable. An increased amount of dark energy should be detectable 
by other cosmological observations. The supernova data alone are not very precise, but 
still homogeneous Friedmann models that fit the data best are the elliptic, not flat, 
models with Q\ 0.9 |52| . Also, the BAO analysis is in the agreement with larger 
amount dark energy. However, it should be noted that one of the elements of the BAO 
analysis is a conversion from the redshift to distance. This conversion is obtained using 
the Friedmann models. Therefore, if the distance-redshift relation is changed then the 
result of the BAO analysis can also be different. This should be kept in mind when 
analyzing the BAO data, especially that as shown in [34j ignoring inhomogeneities also 
leads to additional systematics. 

5. Conclusions 

This paper investigated how the presence of inhomogeneities (small-scale inhomo- 
geneities observed in the Universe) affects the distance to the surface of last scattering, 
and how via the change of the distance the inhomogeneities influence the CMB data. 
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In the standard approach, the distance is calculated within a framework of 
homogeneous Friedmann models. There are several reasons for using this approach. 
First of all, there is an argument used by Weinberg [53j who pointed out that although 
for a single case the distance is modified by the inhomogeneities, but due to photon 
conservation, when averaged over large enough angular scales the overall effect is zero 
(however see [5l] for an argument why Weinberg's reasoning should not apply). Second 
of all, as seen from ((7|) if the amount of voids is the same as amount of overdensities then 
Ad should be zero and the distance should be exactly the same as in the Friedmann 
model. As density fluctuations should vanish after averaging over sufficiently large scales 
it seems reasonable to expect that the Friedmann distance-redshift relation correctly 
describes the reality. 

However, there is a difference between vanishing 3D average and vanishing of the 
average of the density fluctuations over the line of sight. If density fluctuations in the 
Universe are purely random then vanishing 3D average implies vanishing average over 
the line of sight. If density fluctuations are not purely random, i.e. if there is some 
degree of organization, then vanishing 3D average does not necessarily imply vanishing 
density fluctuations over the line of sight (for a discussion see |35]). As the large-scale 
matter distribution in the Universe has a form of the cosmic web, there is some kind 
of organization, and thus the result of ([7]) does not have to be zero, even if 3D average 
of density fluctuations vanishes after averaging over sufficiently large scales. This is the 
case of model Len 2. Also, relation ([7]) does not apply when structures are non-linear. 
The calculations within model SC were carried out within the non-linear regime. In 
addition matter distribution was chosen so that the density fluctuations along the line 
of sight vanish after averaging (i.e. (5) id = 0). In this case, again, the distance was not 
the same as within the Friedmann model. 

Even if the change of the distance is small, like in Len 1 model, still the effect is 
important, as we have already entered the era of precision cosmology - as seen from 
Table [1] and Table [2] the change of the inferred cosmological parameters, even for Len 
1 model, can be of order of a percent or more. If the change of the distance is larger, 
for example as in models Len 2 and SC, then the inferred parameters are significantly 
different. In this case spatially flat models are ruled out and the data favors positively 
curved models with with Q\ 0.8 — 0.9. 

As discussed in Sec. ([3]), within a chosen method, the distance to the last scattering 
surface is relatively uniform with only small variance. The presented, in this paper, 
results show that the change of the mean is important, and the proper handling of 
this effect is essential, otherwise the systematics may be larger than the precision of 
measurements. The only problem, however, is to have a good method for the calculation 
of the distance. All models presented in this paper have their limitations, and should 
rather be treated as toy models, i.e. as examples to show the significance of the problem. 
The results show that the accuracy of the CMB analysis strongly depends on the 
accuracy of the calculation of the distance corrections. Peebles described the importance 
of accuracy and precision in cosmology with the following example a digital scale may 
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read out the weight of an object to many significant figures, in a precise measurement. 
But if the scale is not well calibrated the measurement may not be very accurate [55J. 
In our case, the 'calibration' requires an inhomogeneous framework that will allow us 
to estimate the distance with inhomogeneities taken into account. As the non-linear 
corrections are important this cannot be done within the lensing approximation ((7j) but 
requires a fully inhomogeneous framework. Only then we will be able to accurately 
predict the distance and say: 'the distance is larger by 29.72%' or '36.14% than 
when homogeneity is assumed'. With the increasing precision of CMB experiments, 
in particular Planck, Atacama Cosmology Telescope, and South Pole Telescope this is 
important. Otherwise, what would be the meaning of a measurement with very high 
precision (~ 1% statistical errors) if the accuracy is low (~ 30% systematics)? 
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Appendix A. Construction of model Len 1 

In order to solve ([6]) one needs to know V 2 along the line of sight. Using the Poisson 
equation the gravitational potential is related to density fluctuations. Therefore, a 
method of generating density fluctuations is needed. In the non-linear regime, the 
density fluctuations are not Gaussian (there is no symmetry as — 1 < 5 < oo). However, 
in the non-linear regime the distribution of density fluctuations can be approximated 
with the one-point log-normal PDF |56l 157] 

(ln(l + 5) + al/2f- 



P(5) = -J==e W 



2^ 



1 + 5' 



(A.l) 



where 



a 2 nl = ln[l + all and o\ = JL J dkV(k)W 2 (kR)k 2 , (A.2) 

o 

and V(k) is the matter power spectrum. Therefore, the algorithm for 5 along the line 
of sight is as follows: 

(i) The radius of a structure R is randomly generated from a uniform distribution, 
from to 10 Mpc. 

(ii) Then or is calculated using ( 1A.2|) . 

(iii) An initial value of a density fluctuation 5q is generated from the log-normal 
distribution flA.ll) . 

(iv) The evolution of 5 (at a fixed point) is calculated using the linear approximations 

[58] 

•■ h ■ AnG . , . 

5 + 2-5 = ^ r p5. A.3 

a c z 

(v) Using the Poisson equation, V 2 is calculated and inserted into (J7|) which is solved 
from Xi to x = Xi + x( 2 ^)- 

(vi) Steps (i)-(v) are repeated so (I7j) is solved from x — to Xe- 

A sample of density fluctuations along the line of sight is presented in Figure [TJ 



Appendix B. Construction of model Len 2 

Instead of using directly the log-normal distribution to generate density fluctuations 
along the line of sight an alternative approach is used. In (35] it was shown that 
compensated structure can also have a log-normal PDF. Thus, Len 2 generates density 
fluctuations along the line of sight using a method discussed in [35J. The algorithm is 
as follows: 

(i) First the radius of a void R v is generated from the Gaussian distribution with the 
mean of 12 Mpc and the standard deviation of 2 Mpc. 
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(ii) Density within the void Q v is generated from the Gaussian distribution with the 
mean 0.2Q m and a = 0.27Q m . If a generated in this way fl v is lower than 0.01Q m 
then the generation is repeated. If after 6 times it is still less than 0.01f2 m then Q v 
is generated for a uniform distribution between and 0.01Q m . If fl v > Q m then Q v 
is generated one more time. If after 6 times Q v > Q m then its value is chosen for a 
uniform distribution from 0.85f2 m to Q m . Density contrast is then 5 = Q v /Q m — 1. 

(iii) Density of the surrounding shell Q s is generated from the Gaussian distribution 
with the mean of 1.75Q m and a = 0.7Q m . If Q s < Q m then its value is generated 
again. If after 6 times Q s < Q m then its value is generated for a uniform distribution 
between 1.75f2 m and 1.95f2 m . Density contrast is then 5 = Q s /Q m — 1. 

(iv) The condition, that the structure is compensated implies that the radius of the 
whole structure is 

a - a, ^ 1/3 



R — R r 



I 



The angle at which the light ray enters the structure is generated from a uniform 
distribution between and 0.257T. 



(vi) The evolution of 5 (at a fixed point) was calculated using (1A.3[) . 

(vii) Integral ([7]) is solved from Xi to Xf (where Xi is the comoving coordinate of the 
entry point and Xf the point where the light ray exits the structure). 

(viii) Steps (i)-(vii) are repeated so ([7]) is solved from x = to \ e . 

A sample of density fluctuations along the line of sight is presented in Figure [TJ 
The reason for defining the profile as above is to have a profile that allows for some kind 
of randomness (the size and the density of voids and shells is random), but still because 
we use spheres (which suppose to correspond to real wall that surround real voids) the 
fluctuations are not purely random - the pure randomness is gone once we decided to 
use spheres. Still because of the above construction the 3D density PDF of this model 
is almost log- normal, (see [35j for a discussion). If instead of the Gaussian, the uniform 
PDF is used (steps (i)-(iii)) then the PDF of S is not log- normal PDF. 



Appendix C. Construction of model Len 3 

(i) First the radius of a void R v is generated from a uniform distribution between 
and 30 Mpc. 

(ii) Density contrast of the void 5 V is generated from a uniform distribution between 
—0.5 and —1. 

(iii) Integral ([7]) is solved from Xi to Xf (where Xi ls the comoving coordinate of the 
entry point and Xf the point where the light ray exits the void). 

(iv) Then the length of a filament Lf is generated from a uniform distribution between 
and 5 Mpc. 
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(v) The density contrast of the filament 5f is generated from a uniform distribution 
between 0.5 and 2.5. 

(vi) Integral (J7J) is solved from Xi to Xf = Xi + x(Af)- 

(vii) Steps (i)-(vi) are repeated so (J7J) is solved from x — to Xe- 

A sample of density fluctuations along the line of sight is presented in Figure [TJ 



Appendix D. Construction of the Swiss Cheese model 

The Swiss Cheese model considered here is based on the Lemaitre-Tolman (LT) model 
[491150] . which is the simples generalization of the Friedmann models. The metric of the 
LT model is 



ds 2 



dt 2 - 



R. 



-dr 2 - R 2 (t, r) (dtf 2 + sin 2 tfd<p 2 ) 



(d.i; 



1 + 2E 

where E{r) is an arbitrary function. In the Friedmann limit E — > —k^r 2 (in a general 
LT model the 'curvature index' is position dependent). R(t,r) is the areal distance and 
in the Friedmann limit R — > ra(t) (in a general LT model the scale factor is position 
dependent). The evolution (i.e. the generalized Friedmann equation) is 



R 2 



2M A „ 
2£ + lf + 3 fl ' 



M(r) is an arbitrary function defined by the mass density: 
2M, r 



R?R, r 



From (ID. 21) the age of the Universe, r, is 

H 

r(r) = t - t B {r) = 

o 



dR 



(D.2) 



(D.3) 



(D.4) 



2E + 2M/R + |Ai? 2 

where ts(r) is an arbitrary function and in the Friedmann model ts — >■ const (in a 
general LT model the age of the universe is position dependent). There are three 
arbitrary functions in the LT model ts(r), M(r), and E(r), however only 2 are physically 
independent. For a review of different applications of the LT models see |59l (id. 61j. 
Here we will assume that the model is defined by the assumption that the age of the 
universe is everywhere the same, te = and the function M(r) is 



M = Mn + < 



( M x e 

M 2 exp 



i-2x 







for £ ^ x a , 

for x a ^ £ ^ 3x a 

for 3x a ^ £ ^ 4x a , 
for £ ^ 4x a , 



(D.5) 



where £ = r/kpc, M = (4vrG/3c 2 )p f) £ 3 , p fe = fJ m g, fi m = 0.262, ^ = 71.4 km s" 1 
Mpc" 1 , M x = x~ 3 M 2 e- 1 - 5 , M 2 = 3.29 x 10 10 kpc and -2.91 x 10 11 for regions A and B 
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respectively, 10 3 and 15 x 10 3 , for regions A and B respectively. Region A and 

B are matches alternately at r = 4x a . The above profile was chosen for the following 
reasons: it behaves like a FLRW model for £ ^ x a but with lower density than outside, 
then for x a ^ i ^ 3x a we have a transition region, and a cubic behavior for 3x a ^ £ ^ 4x a 
allows for a smooth matching to the background values. Density fluctuations along the 
line of sight are presented in Figure [TJ It should now be clear why we have selected the 
above set of functions. These functions allow us to obtain the present-day density profile 
as presented in Figure [U This model has thus 2 advantages: 1) it allows to calculate the 
evolution and the distance exactly, without additional approximations, 2) the mean of 
the density fluctuations along the line of sight is zero, i.e. as in Len 1 model and unlike 
in Len 2 and Len 3 models. 

When constructing a Swiss Cheese model, the junction conditions need to be 
satisfied. Here we match the LT inhomogeneity (holes) to the Friedmann background 
(cheese). The LT patches are placed so that their boundaries touch wherever a light ray 
exits one inhomogeneous patch. Thus, the ray immediately enters another LT patch and 
does not spend any time in the Friedmann background. To match the LT patch to the 
Friedmann background across a comoving spherical surface, r = constant, the conditions 
are: 1) the mass inside the junction surface in the LT patch is equal to the mass that 
would be inside that surface in the homogeneous background; 2) the spatial curvature 
at the junction surface is the same in both the LT and Friedmann models, which implies 
that E(r) = —k r 2 and E' = —2k r; 3) the bang time and A must be continuous across 
the junction. As seen, by the construction these conditions are satisfied. 

Finally, the redshift is given by |61j 



Only radial geodesic (i.e. the geodesies that propagate through the origin) will be 
considered. In such a case the shear of light bundle is zero [271 [28] . Thus, the junction 
of null geodesies only requires the matching of one of the components of k a as the other 
one is given from k a k a = 0. 



dr 
di 
dt 

d~z 



A 1 yT+ 2E 
l + z R> 
1 R! 



(D.6) 



